Derivation of symmetry-based pseudo spin models for modulated materials 
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A general concept for the derivation of symmetry-based 
pseudo spin Hamiltonians is described. It systematically 
bridges the gap between the atomistic basis and various 
pseudo spin models presented in literature. It thus allows 
the application of the multitude of analytical and numerical 
procedures derived for the statistical mechanics treatment of 
the latter to the description of structurally modulated crys- 
tals and furnishes a general frame for reviewing these models 
from a common point of view. 



I. INTRODUCTION 

Various approaches for the description of systems 
exhibiting commensurately and/or incommensurately 
structurally modulated phases .were presented in the lit- 
erature™ Aandau- type modelsa'H, the model of Chen and 
WalkerBu for A2BX4-compounds and betaine calcium 
chloride dihydrate (BCCD), the model of Hlinka et al.El 
for BCCD, other microscopic models with continuous lc& 
cal variables like the D1FFOUR (discrete $ 4 ) modelsEl 
and pseudo spin models, like the ANNNI (axial next 
nearest neighbor Ising)l3@, AANMDI (axial antisymmet- 
ric nearest neighbor double Ising)B and D1S (double Ising 
spin)EHrEil models. .Ab-|initio and molecular dynamics 
(MD) investigationsllalij allow to gain insight into pro- 
cesses and interactions on an atomistic level. 

A comparative analysis of the merits and shortcom- 
ings of these models is facilitated and systematized by 
the development of a general procedure for the formu- 
lation of model Hamiltonians for uniaxially structurally 
modulated compounds. Such a procedure is described 
in this paper, ft is based on the discrete arrangement 
of atoms and makes use of the symmetry of the lattice, 
ft leads to a hierarchy of Hamiltonians corresponding 
to different approximations, that describe, among oth- 
ers, the different microscopic models mentioned above. 
These can thus be discussed from a common point of 
view. The models are obtained by identification of the 
relevant variables and specification of the respective cou- 
plings. Their Hamiltonians allow both numerical and 
approximate analytical treatments. One advantage of 
such atomistic symmetry-adapted model Hamiltonians is 
that they allow predictions also on complicated phase di- 
agrams containing long-period structures, on their space 
group symmetries, and on phase transitions which are 
not accessible by ab initio and MD calculations. Previ- 
ously it was considered a weak spot, especially of pseudo 
spin models, that the significance (on an atomistic scale) 



of the model (pseudo spin) variables was not well defined 
and that predictiousjn terms of space group symmetries 
could not be madeO'Ll 



II. FORMULATION IN TERMS OF 
SYMMETRY-ADAPTED LOCAL MODES 

The basis of the procedure is to start from symmetry- 
adapted local modes (SALMs) as a localized basis set for 
the expansion of atomic displacements in properly cho- 
sen crystallographic (sub)cells. This facilitates the incor- 
poration of both the discreteness of the lattice and the 
overall symmetry of the crystal. The relevant SALMs 
can either be determined by first-principles calculations 
or taken from the observed displacements at the phase 
transitions. Generalized local variables defined in terms 
of the amplitudes of these SALMs serve as variables for 
a microscopic Hamiltonian. This procedure stresses the 
atomistic roots of the models, leads to models conform- 
ing to lattice symmetry and facilitates the prediction of 
spontaneous polarizations of modulated phases and of 
the displacements occuring at structural transitions. 

Let the crystal be composed of N unit cells with K 
particles each (of which K' are non equivalent). The 
crystal is described by a 3K A^-dimensional configuration 
vector R = R + <5R, whose entries specify the average po- 
sitions of the atoms. R describes a T- and p-dependent 
reference structure which is invariant under all elements 
Q = {R|t} of the normal phase space group Go- The 
symmetry-breaking displacement (5R is zero in the nor- 
mal phase. In the lower symmetry phases (unmodulated 
or commensurately or incommensurately modulated), it 
can be written as a superposition 



(i) 



of local contributions i5R n having non-zero entries for 
only those atoms which are associated to cell n. 

<5R n is expanded in terms of a local basis set. If the 
set reflects the crystal symmetries, it is sufficient in most 
cases to retain only a few (e.g. one or two) local modes. 
The smallest unit in the crystal with respect to symme- 
try is the asymmetric unit Aq of the space group Go- 
Symmetry imposes no restrictions on the positions of the 
K' atoms associated to this polyhedral subcell (nor need 
they to be known except for a quantitative atomistic de- 
termination of model parameters). Further subcells A„ 
defined by 
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A* - ^A and = {E|n}A 

cover the whole crystal. The B = [Go : T] space 
group operations Q 1 , . . . , Q B are coset representatives in 
a coset decomposition of Go with respect to its sub- 
group T of primitive translations. They have to be cho- 
sen such that the subcells A^,...,A„ form a contigu- 
ous new choice of cell n. The basis set of 3K' subcell 
modes V^, k = 1, . . . , 3K' for the displacements of the 
K' atoms associated with A^ should satisfy the relation 
V„/ K = GY^n: where Q E Go is the space group operation 
transforming A^ into A^,. It is thus entirely determined 
by the choice of modes in the asymmetric unit Aq. 

A complete - for K'B > K over-complete (see below) 
- set of 3K'BN local modes can be constructed as linear 
combinations 

B 

Wl K = E sH ^n (2) 

i=i 

of the modes in subcells i — 1, . . . , B with a B x B ma- 
trix §. The subset {W„ K : k = 1, . . . , 3iT} is called the 
7th set of SALMs for cell n. The transformation matrix 
S should be chosen in accordance with the irreducible 

representations (IRs) of the point group Go such that for 
the construction of every symmetry mode (deformation 
<5R transforming according to an IR of the group G(q) 
of the wave vector) only a minimum of different sets of 
SALMs is needed. Local modes W„ K and describe 
the same displacement patterns in their respective cells, 
i.e. = {E|m - n}W^. All linear combinations of 

displacement vectors in one set of SALMs have the same 
transformation behavior, which is entirely determined by 
S, since, apart from a primitive translation, an arbitrary 
space group operation only permutes the subcell indices. 
Only those sets of SALMs can contribute to the defor- 
mation SK that have the symmetry properties imposed 
by the IR according to which <5R transforms. 

Atoms on special positions (faces, edges and corners of 
the subcells) cause the number SK'BN of local modes to 
be higher than the number 3KN of independent displace- 
ments. Hence, the coefficients of an expansion in terms 
of these SALMs have to satisfy consistency relations. 

The symmetry-breaking displacement £R n in cell n is 
decomposed into contributions <5R^ from the B sets of 
SALMs: 

B B 3K' 

<$Rn=E^ = EE a n K W^. (3) 

7=1 1=1 K =l 

Under a change of external parameters, <5R n will fol- 
low a valley of low local potential (see below) in the 
high-dimensional configuration space corresponding to 
the modes relevant for the phase transitions. In the sim- 
plest and frequent case this valley is not very shallow and 

the contribution <5R„ = J2 K =i a n K W^ K from the Ith set 
of SALMs will lie on a curve (dimension d a = 1). In a 



more general case, <5R^ will lie onarf < 3if'-dimensional 
curved surface E. This generalization would account for 
deviations due to different (small) external fields and al- 
low the treatment of additional fluctuations (about the 
bottom of the valley) described by the additional d a — l 
modes. The special form of the symmetry-breaking dis- 
placement (in the asymmetric unit Aq) connected with 
the transitions is not determined by space group sym- 
metry. It will, in general, be different even for materi- 
als exhibiting the same symmetries of the normal and 
the modulated phases (e.g. A2BX4 compounds). One in- 
troduces local sets of curvilinear coordinates for every 
set of SALMs by means of a coordinate transformation 
ai K = f lK (Q?, Q^ K '), k = 1,...,3K' leading to 

3K' 

K=l 

such that E is parametrized by the first d a coordinates 
(relevant variables per set of SALMs). 

The space group and the direction of the modulation 
determine the subcell generators Q 1 ,.. . , Q B and the ma- 
trix S and thus should be the same for a whole class 
of materials. On the other hand, the displacements of 
atoms in the asymmetric unit and the coordinate trans- 
formations f lK are determined by the local potential and 
thus by the specific material under investigation. This 
allows a separation of effects of crystal symmetry from 
effects characteristic for a special substance. 

The symmetry-adaptation (||) provides a two-way re- 
lationship between the B sets of SALMs and the IRs of 
G(q): 1) In the superposition (||), a symmetry mode <5R 
transforming according to a given IR requires the inclu- 
sion of a few distinct sets of SALMs (e.g. I = 1, . . . , d s ) 
and the exclusion of the others. 2) The included sets, on 
the other hand, may lead, with different spatial modula- 
tions of the respective relevant local variables Q 1 ^, . . ., to 
displacements transforming according to different (but in 
general not all) IRs. 

In the simplest case - one relevant set of SALMs (d s = 
1) and one degree of freedom in the asymmetric unit 
(d a = 1) - expansion (||) simplifies to 

<m = E<^(On\o,...,o) 

n 

and the number of variables per unit cell is reduced from 
3K to 1 (Qi 1 ). 



III. DERIVATION OF MICROSCOPIC 
SYMMETRY-ADAPTED HAMILTONIANS AND 
RELATION TO OTHER MODELS 

The potential is expanded in a power series of the Q^ K . 
The terms containing only variables with given n form 
the local potential $ loc (Q^, . . .,Q B ' 3K ') for cell n. The 
general Hamiltonian of the system is then 
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H = 



E ( p ^) 

I,K 



+ $ int (. 
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where P^ K is the momentum conjugate to <5„ K - Averag- 
ing out all irrelevant generalized coordinates by taking 
the trace over variables I ^ 1 , k ^ 1 yields an effective 
Hamiltonian 



n = tv(pn) = E [(p?) 
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In the case of two relevant sets of SALMs per cell (d s = 
2, I = 1, 2), one gets the two- variable version 
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A form with the superscript 21 replaced by 12 is obtained 
in the case of one relevant set of SALMs per cell (d s = 
I, I = 1) but two relevant generalized coordinates (d a = 
2) per asymmetric subcell. 

The total number of relevant variables per cell is d a -d s . 

The local potential $ and the interaction potential $ 
are effective potentials. Due to the thermal averaging 
over the irrelevant variables, they depend, by their very 
definition, directly on temperature and stresses (which 

both determine the strains) . Usually $ ° c and <i> lnt are 
expanded up to forth and second order respectively and 
only interactions between first, second and third nearest 
neighbors are kept. 

The Hamiltonian (Q) models only the symmetry- 
breaking contribution per cell. The totally symmetric 
IR of Go, which also contributes to the displacements, is 
always present but has been absorbed into the atom po- 
sitions of new reference states, which can be viewed (for 
the whole range of T and p) as (unstable) normal phase 
configurations. 

For modulated crystals with normal phase space group 
Go = Prima exhibiting a pseudoperiodicity along the 
direction of modulation, half cells instead of crystallo- 
graphic unit cells may be chosen. The Hamiltonian (^) 
corresponds to a version of the DIFFOUR model with 
one degree of freedom per lattice site, or to the ANNNI 
model. The case where two local modes (per half cell) 
with different transformation behavior (d s — 2) are in- 
corporated corresponds to Chen and Walker's model or 
to the DIS model. Assuming, in addition, two relevant 
(d a = 2) generalized variables Q 1 ^ and Q 1 ^ for each of the 
two relevant sets of SALMs, one arrives at a model corre- 
sponding to a version of the DIFFOUR model with four 
degrees of freedom per lattice site. If Q 1 ^ and Q 1 ^ are 
assumed to describe only displacements of the betaine 
and|-CaCl2 groups respectively, the model discussed in 
refsJalia for BCCD is obtained. 



IV. INTRODUCTION OF THE PSEUDO-SPIN 
FORMALISM 

Space group operations can be considered to transform 
the generalized variables instead of the modes W „ 
[p = (Ik)). The local potential may have equivalent min- 
ima: if there is a minimum at a certain set of the variables 
Q£ then every space group operation which transforms 
a cell into itself but changes the configuration Q£ gener- 
ates another equivalent minimum. A further simplifica- 
tion may then be introduced by projecting the remaining 
continuous generalized local variables 



Q£ = |Q£|sign Ql=: |Q>£. 



(5) 



onto two-valued pseudo spin variables cr£ representing 
their signs. Because they assume only the values +1 or 
— 1, they are called pseudo spins. They should not be 
confused with usual spins nor are they localized at a def- 
inite position in the cell. A switch of erg from +1 to — 1 
or vice versa corresponds to a transition from one well 
to another and is connected with a collective motion of 
all particles in the unit according to the corresponding 
SALM W„. The motion in a single well as described by 
Q^l occurs on a much shorter timescale and can there- 
fore be averaged out. The interactions (expanded up to 
second order in QPJ are replaced by the averages over 
single wells (taken separately for every |Q£|): 
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with couplings = <f nm • \Qn \ ■ \Qm\- Since the latter 
are defined as thermal averages, they depend inherently 
on external parameters. The procedure yields the follow- 
ing pseudo spin Hamiltonian: 



(6) 



n m pq 



Special systems discussed along these lines were the 
ANNNI model (one rekjaiit mode per unit) or the DIS 
model (two variables)lij'Ll Symmetry considerations 
yield conditions for the dependence of the potential on 
the generalized variables, which have for example been 
taken into account when formulating the AANNDItJ or 
the DISS model. 



V. SYMMETRY OF THE STRUCTURAL 
MODULATION 

Directly below the transition temperature from the 
unmodulated high temperature phase to the modulated 
phases a sinusoidal structural modulation is formed, 
which can be described by a Bloch mode of distinct 
symmetry (symmetry mode). Upon further cooling, it 
becomes more and more rectangular. Thus higher har- 
monics or further symmetry modes need to be consid- 
ered. In good approximation, the relative displacements 
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in one unit (cell or half cell) do not change whereas the 
amplitudes vary from unit to unit. This enters in the 
present approach by using restricted (e.g. one mode) lo- 
cal displacements (SALMs) and modeling the sinusoidal 
modulation and deviations from it by a variation of the 
amplitudes of the SALMs. 

The local displacements are determined by the local 
potentials; hence the relevant SALMs and the respective 
coupling parameters can be determined for example by 
ab initio-type calculations. With the methods of statisti- 
cal mechanics, the resulting Hamiltonians yield the sta- 
ble profiles of the model variables. The superpositions of 
SALMs described by these profiles correspond to a mix- 
ture of a few symmetry modes. Just below the normal 
phase, where the profile is nearly sinusoidal, the primary 
contribution to the structural deformation of the crystal 
transforms according to an IR of G(q). 
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VI. CONCLUSIONS 



Symmetry-adapted local modes defined for crystallo- 
graphic (sub)cells represent a natural lattice theoretical 
basis for the formulation of model Hamiltonians for de- 
scribing modulated systems. Symmetry considerations in 
terms of the group of the modulation vector lead to the 
selection of admissible local mode combinations. By well 
defined thermal averaging, the lattice theoretical systems 
can be projected onto corresponding pscudo spin models. 
The coupling constants of the latter depend, by their very 
definition, on temperature and strains. Explicit calcula- 
tions along these lines would allow, for not too compli- 
cated systems, to go all the way from an atomistic theory 
to explicit statements on phase diagrams, polarizations 
etc. In such a treatment approximations are of course 
necessary and should be made (in the original atomistic 
theory as well as in the consecutive model calculations); 
but one should avoid ad hoc assumptions (fits) as far as 
possible. 
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